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ON 



The phase diagram of ultracold bosons in realistic optical superlattices is addressed via second- 
order cell strong coupling perturbative expansions for the Bose-Hubbard model describing the sys- 
tem. Taking advantage of the cell partition inherent in the complex periodic modulation of a super- 
lattice, this technique allows for the description of the unusual loophole-shaped insulator domains 
that may occur in the phase diagram of the system, unlike the standard perturbative approach. 
Furthermore, comparisons with quantum Monte Carlo simulations show that our approach pro- 
vides quantitatively satisfactory results at a significantly lower computational cost than brute force 
numerical methods. We explicitly consider the phase diagrams for two realistic 3-periodic optical 
superlattices. These show that many insulator domains exhibit an unusual reentrant character, 
which we discuss, and suggest that the quantum phase transition relevant to the loophole-shaped 
insulating domains does not require extreme experimental condition in order to be observed. 

> . 

^ , PACS numbers: 03.75.Lm, 74.81. Fa, 05.30.Jp, 73.43.Nq 

S ■ I. OVERVIEW 

■ 

In a recent experiment 1] , Bose-Einstein condensates (BECs) were loaded in a one-dimensional trapping potential 
created by superimposing two independent one-dimensional homogeneous optical lattices with commensurate lattice 
i constants di — l\a and c?2 = ^2^, <j 6 tf |. This gives rise to an optical trapping potential characterized by a 

periodic modulation involving £ = max(£i,£2) subsequent local minima, i.e. a ^-periodic superlattice. 

In this paper we address the supcrfluid-insulator quantum phase transition that is expected to characterize this 
system, due to the competition between the boson repulsive interaction and the kinetic energy J3( . We recall that the 
latter can be tuned experimentally, being directly related to the strength of the optic potential 4] . The control of this 
parameter actually played a key role in recent experiments where the superfluid-insulator quantum phase transition 
£^ \ was revealed for two and three dimensional arrays of Bose-Einstein condensates 0, . 

In particular, we provide the zero-temperature phase diagram for the Bose-Hubbard (BH) model, originally intro- 
l/"") ' duced for liquid helium in porous media Island nowadays commonly adopted in the description of BEC arrays trapped 
in optic lattices Q and superlattices 0, SH ■ We recall that in the latter case the phase diagram may exhibit unusual 
loophole shaped insulator domain corresponding to fractional fillings. These domains, that cannot be described by 
standard strong coupling perturbative expansions (SCPE), have been shown to exist by exploiting the exact mapping 
between hard core bosons and free spinless fermions on one dimensional systems || , and subsequently investigated by 
means of a qualitative multiple-site mean field approach [Toj . In the case of £ — 2 superlattices, quantitative results 
have been provided in terms of (analytical) first-order cell SCPE and quantum Monte Carlo simulations 0. 

Here, such phase diagram are obtained in terms of second order cell SCPE for the ground state energy of the system, 
that — unlike the standard SCPE — allow a rather satisfactory quantitative description of all of the critical fractional 
fillings characterizing superlattices 0. This is possible basically because in cell SCPE the unperturbed system is not 
a collection of individual lattice sites , but rather the collection of the individual unit cells inherent in the richer 
periodic nature of superlattices. For all of the fractional fillings determining the key features of the zero-temperature 
phase diagram, this choice allows the analytic description of the unperturbed ground state that lies at the basis of a 
perturbative expansion. 

The price for the improved description granted by cell SCPE is the loss of the completely analytical treatment 
allowed by standard SCPE. Indeed, completely analytical results based on cell SCPE are approachable only at first 
perturbative order, for small cell sizes and fillings [9j- On the other hand, cell SCPE provides quite satisfactory results 
with a relatively small numeric effort, consisting basically in the exact diagonalization of the BH Hamiltonian relevant 
to the isolated unit cell of the superlattice. This allows a reliable quantitative description of situations beyond the 
reach of standard SCPE and/or quite demanding for more refined numeric approaches, such as quantum Monte Carlo 
simulations (relatively large fractional fillings for superlattices whose unit cell comprises several sites). 
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The most noteworthy result evidenced by our cell SCPE concerns the critical kinetic to repulsive energy ratio char- 
acterizing the loophole insulating domains. Indeed, according to previous estimates based on a mean-field approach 
, the loophole insulating domains are entered at critical values of the above ratio much smaller than that relevant to 
the usual lobe-shaped domains. Our second-order cell SCPE, which provide a more quantitative estimate than mean- 
field approximations, show that the critical values for the loophole domains are of the same order, and sometimes even 
larger, than those of the usual lobe domains with comparable critical fillings. This suggests that the insulating phase 
corresponding to fractional filling loophole domains does not require extreme experimental conditions. Furthermore, 
second-order cell SCPE confirm the reentrant character of the insulating domains predicted for superlattice by a 
mean- field approach and first-order cell SCPE 

The plan of this paper is the following. In Sec. [n]we recall the superlattice Bose-Hubbard model 0, B @ f° r the 
description of ultracold atoms in optic superlattices and introduce the key ideas of cell SCPE, the details of this 
approach being confined in the appendix for the sake of readability. In Sec. IIIII we illustrate the capabilities of the 
technique by comparing cell SCPE and quantum Monte Carlo results for a I = 2 superlattice. We confirm that cell 
SCPE is satisfactory also when the perturbative parameter is not very small. In particular quite satisfactory results 
are obtained even in the the homogeneous case. Furthermore we provide the zero-temperature phase diagrams for 
some realistic one dimensional superlattices pj. Section iTVl contains our conclusions. 



II. SUPERLATTICE BH MODEL AND CELL SCPE 

Following the tight-binding-like approach of Ref . 0| , ultracold atoms in optical superlattices are commonly described 
by means of a superlattice Bose-Hubbard model (3, la, 111 characterized by a periodic modulation of Hamiltonian 
parameters such as the local energy offsets and the hopping amplitudes between neighbouring lattice sites. In the 
particular case of one dimensional ^-periodic superlattices, a unit cell consists of £ subsequent sites, so that a site can 
be conveniently labeled by two indices, denoting the cell it belongs to and its position within such cell, respectively. 
Hence, the Bose-Hubbard (BH) Hamiltonian for a superlattice comprising C cells each containing £ sites can be written 
in the form 
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joins subsequent cells. Operators a^T c , Ofe, c and n k . c — a k c a k _ c create, annihilate and count bosons at site labeled 
(k, c). As to the Hamiltonian parameters, U > is the boson (repulsive) interaction, \i is the chemical potential, v k is 
the energy offset of the fcth site within the cell, the tk's with k = 1,2, ...,£ — 1 are the intra-cell hopping amplitudes 
and ti accounts for the hopping of bosons across neighbouring cells. Of course on an infinite ^-periodic structure the 
£ possible choices for the boundaries of the cells are in principle equivalent. However, in view of our strong coupling 
perturbative expansion, the most reasonable choice is such that the inter-cell hopping amplitudes are the smallest. 

As we recalled above, the one-dimensional ^-periodic superlattice Bose-Hubbard (BH) exhibits an insulator- 
superfluid quantum phase transition at any rational filling of the form f = L/t, I € N 0. These phase transitions 
are commonly put into evidence in the zero-temperature n/U-t/U phase diagram, where t is a scaling factor for the 
superlattice hopping amplitude related to the strength of the optic potential (that is tfc = t ■ Tfe). More to the point, 
the boundaries of the insulating domain relevant to the critical filling / are obtained comparing the ground state 
energies of the systems containing N = fM, N + 1 and N — 1 bosons, where M is the number of sites in the lattice, 
which should be large enough to provide a satisfactory approximation of the thermodynamic limit. Since the size 
of the Hilbert space for N bosons on M sites is d(N,M) — ( N+ n~ 1 )> exac t diagonalization is beyond the reach of 
common computational capabilities already for relatively small systems at low fillings, and a numeric estimate of the 
ground state energies under concern requires complex algorithms such as quantum Monte Carlo [T2I Il3| or density 
matrix renormalization group simulations |l4l Il5j . 
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In Ref. |ll| j satisfactory approximations of the ground state energies required for the description of the insulator 
domains are obtained for homogeneous lattices in terms of analytic third order expansions in the ratio of the hopping 
amplitude and the boson interaction. This approach can be extended to structures featuring inhomogeneous hopping 
amplitudes, including superlattices |16| . but in general for integer critical fillings only. Indeed, in the case of fractional 
critical fillings, an explicit form for the structure of the unperturbed ground state is might not be available for the 
boson populations involved in the determination of the insulator domain boundaries. This typically happens when 
some of the local potentials Vk within the unit cell are equal, and loophole insulator domains appear in the phase 
diagram Consider for instance the extreme case where there is no modulation in the local potentials (i.e., without 
loss of generality, Vk — 0, Vfc) and the periodic pattern of the superlattice is due to the hopping terms alone, which 
do not appear in the unperturbed Hamiltonian of standard SCPE. The only fractional fillings for which the form 
of the unperturbed ground state can be easily provided are those differing by a single boson from an integer filling, 
i.e. those involved in the determination of the integer-filling insulator domains [ill Il6j . This is basically due to 
the fact that the degeneracy of the lower eigenspace of the unperturbed Hamiltonian equals the lattice size, so that 
finding the unperturbed ground state reduces to a trivial single-particle problem. For any other fractional filling, the 
degeneracy of the lower unperturbed eigenspace is much larger than the lattice size, and the ground state cannot be 
found analytically. 

This problem can be circumvented adopting a perturbative approach that takes advantage of the periodic modulation 
of the Hamiltonian parameters characterizing a superlattice BH model, i.e. whose perturbative parameter is the inter- 
ccll hopping amplitude [9| . That is to say, in the cell SCPE the unperturbed Hamiltonian is the sum of the independent 



cell Hamiltonians 



while the perturbative term is that defined by Eq. ■ We mention that a similar approach was adopted in the study 
of the Hubbard and Heisenberg [l^ models in the special case of a dimerized chain (i.e. a i = 2 superlattice). 
In some sense, cell SCPE is standard SCPE applied to the homogeneous lattice made up by unit cells instead of 
individual lattice sites. This choice makes it easy the determination of the analytic form of the unperturbed ground 
state for the total populations N — CL, N = CL ± 1 relevant to the fractional critical filling / = L/t. Indeed, a 
generic eigenstate and the corresponding eigenvalue of the unperturbed Hamiltonian Q have the form 



C 

|k,N)) =(g)|fc c ,iV c ) c , S(k,N) = 53<« (5) 

c— 1 c 

where k = (fc 1; k 2 , ■ ■ ■ , fee), N = (JVi, N 2 , . . . , Nc), \k, N) denotes the fcth eigenstate of the cell Hamiltonian for a 
population of N bosons, and is the relevant eigenvalue. In the following we assume that these states are sorted 
according to their energy, so that the label k = 1 refers to the ground state. Thus, in particular, the ground state of 
Hq relevant to the critical fractional filling / — Ljl (i.e. to the total population TV = CL) and the relevant energy 
have the form 

|*cz» = (g)|l,£>c, £ ( cl=CEt (6) 

c=l 

Taking into account the translational symmetry of cell lattice, the defect ground states \H\ of Hq, corresponding to 
particle and hole total populations N = CL ± 1 , have the form |19| 
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According to standard perturbative theory, the subsequent terms in the expansion of the ground state energies relevant 
to the total population TV 



= (9) 
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where the prime signals that states with the same energy as the ground state must be excluded (recall that in the 
"defect" cases, N — CN ± 1, the ground state of H is C-fold degenerate [Tsjl- 

As it is discussed in some detail in the Appendix, the perturbative corrections are given in terms of matrix 
elements (on the unperturbed eigenbasis) of the single boson operators at the cell boundaries: 



Generally, such matrix elements must be evaluated numerically. However, the necessary computational effort is easily 
affordable even for cells comprising several sites, and not only at the lowest critical fillings. Indeed the second order 
cell SCPE result for the domain boundary relevant to the critical filling Ljl of a ^-periodic superlattice requires the 
complete diagonalization of a £-site cell BH Hamiltonian © for the five total populations ranging from N = L — 2 
to N = L + 2. For instance, the most demanding operation to be faced when investigating the critical filling 11/5 
for a period 5 superlattice is the complete diagonalization of a matrix whose size is d(13, 5) = 2380, which requires a 
negligible amount of time on a reasonably up-to-date personal computer. 

Note that, in order to obtain more refined results for the same system, one can resort to more complex numerical 
methods, such as quantum Monte Carlo simulations. However, the necessarily finite realization of the superlattice 
addressed by the simulation must comprise a sufficiently large number of cells to provide a good approximation of 
the thermodynamic limit. Since one expects that such number is of the order of the number of sites in a reasonably 
good finite-size realization of the one-dimensional homogeneous lattice, it is clear that even superlattices with a small 
periodicity are computationally quite demanding. 

The explicit form of Eqs. © and (|10fl in terms of matrix elements in Eq. l|llf) and the details of their derivation 
is given in the Appendix. Comparison with finite size exact results confirmed the correctness of our perturbative 
expansions. 



In this section the cell SCPE technique is adopted in the study of the insulator domain borders for several one 
dimensional superlattices. These borders are standardly obtained as fi± — ±[£cl±i — £cl) and plotted in the n/U-t/U 
plane, where we recall that t is a scaling factor for the superlattice hopping amplitudes (i.e. tk = t • Tfe). For a given 
value of t we obtain the ground state energies ScL+a, with a — 0, ±1 from the second-order SCPE in Eq. (JSJ). Notice 
that the perturbative order refers to the parameter ti = tri, and that the unperturbed Hamiltonian, Eq. (0) depends 
on the overall hopping factor t as well. This means that the perturbative estimates of the ground state energies and 
the insulator domain boundaries thereby are not second order functions of t. 

The simplest one-dimensional superlattice is the dimerized chain, corresponding to I = 2 [^. Il7lll8| . As we mention, 
when <i = t-2. the half-integer-filling insulating domains exhibit an unusual loophole shape. A qualitative description 
of these domains can be obtained resorting to a multiple-site mean-field approximation [l(j , whereas second order cell 
SCPE provides quite satisfactory quantitative results. This can be appreciated in Fig. ^ where the domain borders 
of the half-filling insulating domain as obtained from cell SCPE and quantum Monte Carlo simulations are compared. 
One might expect the perturbative technique to be valid in the limit ti -C t\. Note however that inset B of Fig. 
shows that it gives satisfactory results even for relatively large t^jix- 

Actually, cell SCPE can be used to address the homogeneous lattice, setting t% = t\ = t. In this situation, the 
method correctly predicts the absence of fractional-filling loophole insulator domains, whereas the results for the 
integer filling domains are comparable to those provided by third-order standard SCPE , as it is shown in Fig. |3 
Of course, the cell SCPE description of the homogeneous lattice can be performed considering larger "unit cells". In 
general, the thermodynamic convergence of the result does not seem to be very fast with increasing cell size. However, 
we notice that cells as small as four sites provide results exhibiting the expected reentrant character of the insulating 
domain (see inset of Fig- - Furthermore, the comparison between the I — 2 first and second order cell SCPE results 
for the homogeneous lattice suggests that an improvement is more likely produced by an increase in the perturbative 
order rather than in the cell size. In this respect, one might observe that cell SCPE are not really necessary in the 
homogeneous case. Still, as we discuss in the previous section, they are called for when dealing with non trivial 
periodic lattices. 

As a further illustration of the cell SCPE technique, we address two superlattices with the same periodicity as the 
one realized experimentally by the authors of Ref.JlJ, i.e. I = 3. We provide their zero-temperature phase diagrams 
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III. RESULTS 
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FIG. 1: Half-filling insulator domain boundaries for a £ — 2 superlattice BH model with equal local potentials, v 1 = t>2 (the 
corresponding optic potential is sketched in inset A). Main plot: ti — 0.2ti — 0.2t. Inset ti = 0.6ti = 0.6t. Solid lines: 
second order cell SCPE; Dashed lines: first order (analytic) cell SCPE 9]; Errorbars: quantum Monte Carlo simulations on a 
superlattice comprising C = 50 cells. 
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FIG. 2: First Mott lobe for the homogeneous ID lattice according to: first-order, I — 2 cell SCPE ( black dashed thin lines); 
second-order, £ = 2 cell SCPE (black solid thick lines), standard third-order SCPE |7Hl (black dotted thin lines); density matrix 
renormalization group simulations fl^l (best known result, gray solid thick lines). The black solid thin line in the inset is the 
second order £ = 4 cell SCPE result. Note that it exhibits the reentrant character typical of the best result, unlike the standard 
third-order SCPE curves. 



considering insulating domains up to critical filling 7/3. In this respect we mention that, in this particular case, cell 
SCPE would provide results for much larger fillings at a negligible computational cost. Conversely, the same results 
would be very demanding for more refined numerical methods, such as QMC or density matrix renormalization group 
simulations. 

The superlattice considered in Fig. |3 (sketched in the inset) is characterized by even hopping amplitudes (t\ =t% = 
<3 = t) and local potentials v\ — v% — 0.5U and v-i — 0. As expected 0, the critical fillings of the form / = (2 + 3fc)/3 
correspond to loophole insulator domains, while the remaining fillings yield usual lobe shaped domains. Figure 01 
refers to a I = 3 superlattice with even local potentials (v/- = 0) hopping amplitudes tx = = t and £3 = 0.3i. Again, 
cell SCPE provides the expected shape for the insulating domains, only the ones relevant to integer filling exhibiting 
the usual lobe shape . 
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FIG. 3: Zero-temperature phase diagram for a £ — 3 superlattice with ii = — tz — t, Vi = V3 = 0.5U and V2 — (the 
corresponding optic potential is sketched in the inset). The critical fillings characterizing each insulator domain are reported. 
The vertical dashed lines signal some values of the chemical potential where the reentrant behaviour of the insulator domains 
is evident (see the discussion in the text). 



The first striking feature emerging from the phase diagrams in Figs. [3] and 0] is the size of the loophole insulator 
domains as compared to the usual Mott lobes. Note indeed that the mean- field results for the same systems predicts 
much smaller loophole domains 0, llflj . Since SCPE provide a better quantitative estimate than mean-field, the 
quantum phase transition corresponding to loophole domains is expected to occur at not exceedingly small values 
of t/U. Note furthermore that the reentrant character put into evidence in Ref. [l4| for the first Mott lobe of the 
homogeneous one dimensional lattice is pervasive in the phase diagrams displayed in Figs. and ^ That is to say, 
values of the chemical potential \i exist where an increase in the hopping amplitude (i.e. in the kinetic energy) turns 
the system from superfluid to insulating. For instance in Fig. lathis happens at fi = 1.3, fi — 0.2 and even twice 
at fi — —0.2. This can be explained taking into account that an increase in t/U at fixed fi/U can be achieved only 
by increasing the total boson population. A maybe more experimentally relevant reentrant feature is that inherent 
in the loophole shape of some insulating domains. Note indeed that the width of an insulating domain, i.e. the 
chemical potential gap, is related to the excitation gap of the system at critical filling. An estimate of the latter, 
possibly achieved as described in Ref. [20|. would hence allow to discriminate loophole and lobe domains. Indeed the 
excitation gap is expected to decrease monotonically for the latter, while it features a maximum for the former. 



IV. DISCUSSION 



In this paper we introduce the cell SCPE technique for the study of one-dimensional superlattice BH models. This 
technique provides a quantitative estimate of the zero-temperature phase diagram of the model at a low computational 
cost. Indeed we recall that the data presented in Figs. HUH are obtained in a few minutes of CPU time on a standard 
personal computer. Furthermore, cell SCPE provide a description of the insulator domain boundaries also for those 
critical fillings where standard perturbative approaches are ruled out, i.e. those corresponding to loophole domains 
9]. 

In Section IIIII The capabilities of our approach are illustrated by the zero-temperature phase diagrams of some 
experimentally realistic one-dimensional superlattices 0] . We once again emphasize the fact that the loophole domains 
are entered at relatively large values of the parameter t/U driving the quantum phase transition, which suggests that 
the observation of the latter does not require extreme experimental conditions. 

With a reasonable effort, the technique — described in some detail in the Appendix — can be adapted to more 
complex superlattices, such as the trimerized kagome lattice whose experimental realization in terms of optic poten- 
tials is proposed in Ref. 0. We emphasize that on such superlattices the computational demand of more refined 
numeric approaches, such as density matrix renormalization group and Quantum Monte Carlo simulations, can be 
overwhelming. Conversely, the computational cost of second-order cell SCPE does not increase significantly for struc- 
tures such as the kagome lattice, and the only complication arises from the more complex topology of the cell lattice. 
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This simply means that some more terms like those described in the Appendix have to be evaluated. The study of 
the superfluid-insulator quantum phase transition for the BH model on more complex superlattices, where the exact 
mapping with free spinless fermions does not apply, is the subject of a forthcoming paper |2l| . 
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APPENDIX: DETAILS OF THE CELL SCPE APPROACH 



In this section we give a sketch of the derivation of the explicit form of the perturbative corrections in Eqs. © and 
(|10|l in terms of matrix elements of the boson lattice operators at the cell boundaries, Eq. I jllj l. 

Since the perturbative term V, Eq. @, is the sum of terms connecting neighbouring cells, the brackets in Eqs. © 
and IjlUfl are the sum of terms like the following 

((k',N'| a t jaij+1 |k,N)) = j{k'v N j\a\,j\ k i' N j)i ■ j+i<^+i,^+i|aw-H|Ai+i.Jy,-+i>i+i 

• J] c {k' c ,K\I\k c ,N c ) c (A.l) 

and its complex conjugate, where we recall that |k, N)), introduced in Eq. denotes a generic eigenstate of the 
unperturbed Hamiltonian PJl. Putting all of the above matrix elements in a stack, the previous product can be 
represented as 
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where the central column pertains to the operator term, while the lateral columns represent the states involved in 
the product. Each level of the stack refers to a given cell of the superlattice, and contains the corresponding term of 
the product in Eq. I|A.1(I . Of course, the product in the rows containing an identity operator vanishes unless k' = k, 
N' = N, in which case it equals one. Likewise the rows containing the boson operators vanish unless N' = N ± 1, 
but in general there is no constraint on the energy labels k and k'. Note that we dropped the cell labels in the boson 
operators appearing in the stack. Henceforth we drop the cell labels unless strictly necessary. 

We are now ready to evaluate the first order terms of the energy expansions. When the total population in the 
superlattice is J\f — CL, where C is the total number of cells and L is an integer, the pictorial representation of the 
terms in Eq. @ is 
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According to the previous discussion, all of these term vanish. Hence E CL — 0. 
When Af = CL + 1 it is clear that the only non- vanishing terms are of the form 
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and its complex conjugate. More precisely, two of such terms are contributed by each of the C cells. Taking into 
account the normalization factor of the symmetrized ground state \^cl+i)), Eq. (J2J), one gets 

4l+i = (l,i + l|a]|l, J L)(l,L|a 1 |l,L + l) + <l,L|^|l ! L+l)(l,L+l|4|l ! L) 

= 2-(l,i + l|a]|l,i)(l,i|ai|l,i + l) (A.5) 

where in the second equation we take advantage of the fact that the operators and the states are real. Likewise 

4 1 2_ 1 = 2(1, L\a\\l, L — L — l|oi|l, L) (A.6) 

Let us now address the second order term in the cell SCPE, Eq. I|10|) . Each term in the corresponding sum contains 
a product two matrix elements of the the form (|A.lfl . However, most of these product vanish. In order to single out 
the non-vanishing products we adopt a trivial generalization of the stack representation introduced above. The stack 
representing the product under concern is simply obtained placing two stacks of the form (| A. 2|> side by side. According 
to Eq. I|1U|) . the innermost columns of such six-column stack refer to the same unperturbed state. The outermost 
column refer once again to the ground state, or more precisely to one of the possibly degenerate unperturbed states 
appearing in its expression, Eqs. and (|7|l. 

When Af = CL, the unperturbed ground state is a single state of the Fock space, Eq. ©. The corresponding stack 
representation for the second order correction to the ground state energy is 
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By inserting the matrix elements corresponding to stack (|A.7|I in Eq. i|l(J|) one gets 



(2) J ^+ lJ ^ 1 \(j,L + l\a]\l,L)(h,L-l\ ai \l,L) + (j,L + l\a\\l,L){h,L-l\a e \l,L) 
£cL = C f- 1 h 2E[-E^-E^ 



(A. 



where the first term in the numerator corresponds to the process illustrated in stack l|A.7(l , the second to its complex 
conjugate and Dl — d(L,£) is the size of the Fock space for L bosons on £ sites. The overall factor C takes into 
account the fact that the operator V is obtained as the sum of C terms. 

Let us now consider the case where the total number of bosons is j\f = CL. The unperturbed ground state is 
now given by the superposition of Fock states in Eq. 0. This gives rise to five different contributions whose stack 
representations are discussed below. 

First of all we notice that there is a stack giving a contribution of the form (| A. 8|) . yet with a different overall factor. 
This corresponds to 
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producing a contribution of the form 

D J^ D ^ \(j,L + l\al\l,L)(h,L-l\ ai \l,L) \ (j,L + l\a\\l,L)(h,L-l\a e \l,L) 

4 2 2 + m = (c - 2) x: E 
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2E{ - E, 



L+l _ E L-1 



(A.10) 



where the factor C — 2 takes into account the number of terms in V giving rise to the processes described by the stack 
Another contribution corresponds to the stack 















1,L 


I 


1,L 


1,L 


I 


1,L 


I, L+l 


I 


J,L + 1 


J, L + l 


a] 


1,L 


1,L 


oj 


h,L-l 


h,L-l 


ax 


1,L 


1,L 


ai 


f,L + l 


f, L + 1 


I 


1,L+1 


1,L 


I 


1,L 


1,L 


I 


1,L 















(A.ll) 



that vanishes unless j = j' = 1. This produces a contribution of the form 

(2) ^ (l,L + l| a ||l,£)(l,£| a t|fe,L-l)(/i,£-l| ai |l,L)(l,L| ai |l,L + l) 

t CL+\.b 2.^ 0P L rpL+l pL-1 

h=l ZK l ~ h l ~ h h 

l,L + l\a\\l,L)(l,L\a\\h,L- l){h,L- l\a e \l, L)(l, L\a e \l, L + 1) 



E 



h=l 
h=l 



2E{ ~ E^ +1 -Ef- 1 



%L + \\a\\\,L)(\,L\a\\h,L- l)(h,L- l\a 1 \l,L)(l,L\a 1 \l,L + 1) 
2E{ - E[ +1 - E^' 1 



(A.12) 



where the second term is obtained by complex conjugation. 
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The third stack to be considered is 















1,L 


I 


l,L 


l,L 


I 


1,L 


1,L 


ag 


J,L + 1 


J,L + 1 


a\ 


l,L 


1,L + 1 


A 


h,L 


h,L 


ai 


1,L + 1 


1,L 


i 


l,L 


1,L 


I 


1,L 















(A.13) 



where j = h = 1 must be discarded since in this case the innermost vectors belong to the lower eigenspace of Hq. The 
process described by the stack (jA. 131) and its complex conjugate produce the contribution 



r(2) 

'CL+l,c 



Dl+i D l 

E 12 a i> h 

3 = 1 h=l 



(j,L + l\al\l,L)(h,L\ ai \l,L + l) + (j,L + l\4\l,L)(h,L\at\l,L) 



E{ + E{ +1 - E 



J 3 a h 



(A.14) 



where <jjh equals zero if j = h = 1 and one otherwise. 
The next stack we consider is 















1,L 


I 


l,L 


1,L 


I 


1,L 


1,L + 1 


a\ 


h',L 


h',L 


I 


1,L 


1,L 


ai 


J, L + 1 


3-L + l 


a\ 


1,L 


1,L 


I 


h,L 


h : L 


a x 


1,L+1 


l,L 


I 


1,L 


1,L 


I 


l,L 















(A.15) 



which produces a non- vanishing contribution only if h = h! = 1. Furthermore j — 1 must be discarded, since in this 
case the innermost vectors belong to the lower eigenspace of Hq- Therefore, one gets 



Dl + i 

°CL+IA " / j 

+ E 



(l,L + l|q||l,L)(l,L|q 1 |j,L+l)(j,L + l| a ||l,L)(l,L| ai |l,L + l) 

E L+1 _ E L+1 

{l,L + l\a\\l,L)(l,L\a e \j,L + l){j,L + l\a{\l,L)(l,L\a e \l,L+l) 

E^ +1 - E L+i 



3=2 -1 ~3 

D ^ (l,£+l|at|l,£)(l,£| ai |j,£ + l)Q-,L + l| a ||l,L)(l,L| ai |l,£ + l) 



j=2 



E L+1 _ E L + 1 



Finally, the stack 















l,L 


I 


1,L 


1,L 


I 


1,L 


l,L 


ae 


j,L + 2 


j,L + 2 


a\ 


1,L + 1 


1,L+1 


a\ 


h,L-l 


h,L-l 


ai 


1,L 


l,L 


I 


1,L 


1,L 


I 


1,L 















(A.17) 



produces the contribution 

DL ^ DL -A(j,L + 2\al\l,L + l)(h,L- l|oi|l,i) 2 + (j, L + 2\a\\l, L) (h, L - l\a £ \l, L) 



c(2) 
°CL+l,e 



EE 

j=l h=l 



E[ + E[ +1 - Ef +2 - El 



L-l 



(A.18) 
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Putting all together, the second-order correction to the ground-state energy, Eq. (|l(Jfl . is 



,(2) _ -(2) 
CL+1 ~ °CL+l,a 



+ £ (2) + £ (2) +£ (2) +£ 

' °CL + l.b ' °CL+l,c ' °CL+l.d ' °i 



CL + l,e 



(A.19) 
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